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Abstract 

We analyzed the performance of a perturbation theory for nonlinear cosmological dynamics, 
based on the Lagrangian description of hydrodynamics. In our previous paper, we solved hydrody- 
namic equations for a self-gravitating fluid with pressure, given by a polytropic equation of state, 



using a perturbation method. Then we obtained the first-order solutions in generic background 



universes and the second-order solutions for a wider range of polytrope exponents. Using these 
results, we describe density fields with scale-free spectrum, SCDM, and LCDM models. Then 
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we analyze cross-correlation coefficient of the density field between N-body simulation and La 

o 

grangian linear perturbation theory, and the probability distribution of density fluctuation. From 

m 

our analyses, for scale-free spectrum models, the case of the polytrope exponent 5/3 shows better 
performance than the Zel'dovich approximation and the truncated Zel'dovich approximation in a 
quasi-nonlinear regime. On the other hand, for SCDM and LCDM models, improvement by the 

& '• 

effect of the velocity dispersion was small. 
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I. INTRODUCTION 



The Lagrangian approximation for structure formation in cosmological fluids provides 
a relatively accurate model even in quasi-linear regime, where density fluctuation becomes 
unity. The Zel'dovich approximation (hereafter ZA) , a linear Lagrangian approx- 

imation for dust fluid, describes the evolution of density fluctuation better than Eulerian 
approximation P, (|. Although ZA gives an accurate description until quasi-linear regime 
develops, ZA cannot describe the model after the formation of caustics. In ZA, even after 
the formation of caustics, the fluid elements keeps moving in the direction set up by the 
initial condition. ZA cannot describe compact and high density structures such as pan- 
cakes, skeletons, or clumps, while N-body simulation shows the presence of clumps with a 
very wide range in mass at any given time |7(. In general, after the formation of caustics, 
hydrodynamical description becomes invalid. 

In order to proceed with a hydrodynamical description without the formation of caustics, 
the qualitative pressure gradient H and thermal velocity scatter in a collisionless 

medium had been discussed. After that, the 'adhesion approximation' ll| was proposed 
from the consideration of nonlinear wave equations like Burgers' equation. In the 'adhe- 
sion approximation', an artificial viscosity term is added to ZA. As another modification, 
a Gaussian cutoff is applied to the initial power spectrum of density fluctuation and then 
evolves with the ZA. This modified approximation is called the 'Truncated Zel'dovich ap- 
proximation' (hereafter TZA) These modified approximations successfully avoid 
the formation of caustics and describe the evolution for long time. However, the physical 
origin of the modification is not clarified. 

We reconsider the basic, fundamental equation for the motion of matter. The collisionless 
Boltzmann equation describes the motion of matter in phase space. The basic equa- 
tions of hydrodynamics are obtained by integrating the collisionless Boltzmann equation 
over velocity space. In past approximations, such as ZA and its modified version, velocity 
dispersion was ignored. Buchert and Dominguez argued that the effect of velocity dis- 
persion become important beyond the caustics. They also argued that models for large-scale 
structure should rather be constructed for a flow which describes the average motion of a 
multi-stream system. Then they showed that when the velocity dispersion is still small and 
can be considered isotropic, that gives effective 'pressure' or viscosity terms. Furthermore, 
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they argued the relation between mass density p and pressure P, i.e. an 'equation of state'. 
If the relation between the density of matter and pressure seems barotropic, the equation 
of state should take the form P oc p 5 / 3 . Buchert et al. showed how the viscosity term is 
generated by the effective pressure of a fluid under the assumption that the peculiar accel- 
eration is parallel to the peculiar velocity; Dominguez 0, clarified that a hydrodynamic 
formulation is obtained via a spatial coarse graining in a many-body gravitating system, 
and the viscosity term in the 'adhesion approximation' can be derived by the expansion of 
coarse-grained equations. 

Dominguez |Lj| also reported on a study of the spatially coarse-grained velocity dispersion 
in cosmological N-body simulations. The analysis showed that polytrope exponent becomes 
7 ~ 5/3 in a quasi-nonlinear regime, and 7 ~ 2 in a strongly nonlinear regime. Dominguez 
and Melott 20] discussed polytrope exponents of velocity dispersion in N-body simulations. 
According to their results, the exponents depend on a model of initial density fluctuation. 

From these points, the extension of Lagrangian perturbation theory to cosmological flu- 
ids with pressure has been considered. Actually, Adler and Buchert [2^ have formulated 
the Lagrangian perturbation theory for a barotropic fluid. Morita and Tatekawa [3| and 
Tatekawa et al. [23] solved the Lagrangian perturbation equations for a polytropic fluid up 
to second order for cases where the equations are solved easily. Hereafter, we call this model 
the 'pressure model'. 

In this paper, we analyze the density field which is described by the Lagrangian ap- 
proximations; ZA, TZA, and first-order pressure model solutions. We calculate the cross- 
correlation function of density field between the Lagrangian approximation and N-body 
simulation. Furthermore we analyze the probability distribution function of density fluctua- 
tion for confirmation. From these results, we determine a polytrope exponent in the equation 
of state. From our analyses of the cross-correlation coefficient and probability distribution 
function of density fluctuation, we find that the value of the polytrope exponent seems to be 
5/3 for quasi-nonlinear evolution, as Buchert and Dominguez argued ^|. However, for the 
determination of a proportion fixed number in equation of state, we must consider further 
physical processes or carry out high-resolutional N-body simulation. 

This paper is organized as follows. In Sec. m we present Lagrangian perturbative so- 
lutions: In Sec. Ill Al we show the first-order solution of pressure model in the Einstein-de 
Sitter background. For comparison, in Sec. Ill Bl and III C\ we show the solution of ZA and 
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the procedure of TZA. 

In Sec. IHIl we compare the density field between the Lagrangian approximations and 
N-body simulation. In Sec 1111 Al we calculate the cross-correlation coefficient of the density 
field. Though it seems that we can reach a conclusion in this analysis, it is insufficient. 
Therefore in Sec. 1111 Bl we analyze the probability distribution function of density fluctua- 
tion. In Sec. IIVI we discuss our results and state conclusions. 



II. LAGRANGIAN APPROXIMATIONS IN GRAVITATIONAL INSTABILITY 
THEORY 

A. First-order solutions of the pressure model 

In this section, we present perturbative solutions in the Lagrangian description. The 
matter model we consider is a self-gravitating fluid with mass density p and 'pressure' P, 



which is given by the presence of velocity dispersion. T 



re 'pressure' we adopt here is the 



same as was introduced by Buchert and Dominguez ^|, i.e. the diagonal component of 
the velocity dispersion tensor when the velocity dispersion is assumed to be isotropic in the 
Jeans equation jl4| . In Lagrangian hydrodynamics, the comoving coordinates x of the fluid 
elements are represented in terms of Lagrangian coordinates q as 



x 



s(q,t), (1) 



where q are defined as initial values of x, and s denotes the Lagrangian displacement vector 
due to the presence of inhomogeneities. From the Jacobian of the coordinate transformation 
from x to q, J = det(dxi/dqj) = det(5^- + dsi/dqj), the mass density is given exactly as 

P = P^ 1 . (2) 

We decompose s into the longitudinal and the transverse modes as s = V q S + S T 
with V q • 5 T = 0. In this paper, we show an explicit form of perturbative solutions in only 
Einstein-de Sitter universe. For a generic background universe, we obtained the perturbation 



solutions in our previous paper [23 1. 



The transverse modes can be solved easily. The first-order solutions become as follows: 

5 T oc a , a- 1 ' 2 . (3) 



Because the solutions do not depend on 'pressure', the solutions of this mode in ZA become 
the same in form. The transverse modes do not have a growing solution in a first-order 
approximation. 

For the longitudinal modes, we carry out Fourier transformation with respect to the La- 
grangian coordinates q. If we assume a polytropic equation of state P = up 1 with a constant 
k and a polytrope exponent 7, we can write an explicit form of first-order perturbative solu- 
tions. In the Einstein-de Sitter background, the solutions are written in a relatively simple 
form. They are, for 7 7^ 4/3, 

S(K, a) cx a"V4 J±5/{8 _^ (^^^ a^^j , (4) 

where J v denotes the Bessel function of order u, and for 7 = 4/3, 

S{K, a) oc a -i/4±V25/i6-c 2 | J K|V2C 1 ? (5) 

where C\ = 47rGp b (a in ) a^/3 and C 2 = «7/>b( a m) 7_1 ■ Pb an d K mean background 

mass density and Lagrangian wavenumber, respectively. a in means scale factor when an 
initial condition was given. 

Here we notice the relation between the behaviors of the above solutions and the Jeans 
wavenumber, which is defined as 




4:7iGp h a 2 V 
dP/dp(p h )J 



The Jeans wavenumber, which gives a criterion for whether a density perturbation with a 
wavenumber will grow or decay with oscillation, depends on time in general. If the polytropic 
equation of state P = up 1 is assumed, 

K 3 = ^« (37 - 4)/2 ■ (6) 

Equation (jUJ) implies that, if 7 < 4/3, Kj will be decreased and density perturbations with 
any scale will decay and disappear in evolution, and if 7 > 4/3, all density perturbations will 
grow to collapse. We rewrite the first-order solution Equation (J3J with the Jeans wavenum- 
ber: 

S(K, a) oc J ±5/{8 ^ (^|^) ■ (7) 
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B. Zel'dovich approximation 

Z A was obtained as first-order solutions with dust fluid in Lagrangian description [lj . The 
solutions are obtained from the solutions of the pressure model as a limit of weak pressure. 
For example, in E-dS model, when we take a limit k — > in Equations (j3J) and the 
solutions converge to those of ZA: 

S(K, a) oc a, a~ 3 / 2 . (8) 

ZA is known as perturbative solutions which describe the structure well in quasi-linear 
scale. However if the caustics appear, the solutions no longer have physical meaning. 



C. Truncated Zel'dovich approximation 

During evolution, the small scale structure contracts and forms caustics. Therefore if 
we introduce some cutoff in the small scale, we will be able to avoid the formation of 
caustics HQ. In TZA, for the avowee of caustics, we deduce a Gauss.au cutoff to 
the initial density spectrum as follows: 

V(k, t in ) — > V(k, t in ) exp (~k 2 /k NL ) , (9) 

where k^L means 'nonlinear wavenumber', defined by 

„ C^NL 

l = a(t) 2 V(k,t in )dk. (10) 

Jko 

The 'nonlinear wavenumber' depends on the scale factor. The relation between the Jeans 
wavenumber Kj and the nonlinear wavenumber k^L will be discussed in Sec. II V I 



III. COMPARISON BETWEEN N-BODY SIMULATION AND LAGRANGIAN 
APPROXIMATIONS 

In this section, we show a comparison between N-body simulation and Lagrangian ap- 
proximations with two statistical quantities. In our previous paper 2^ , we showed that the 
effect of second order perturbation was still small just before shell-crossing. Therefore we 
consider only first order perturbation. 
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We analyze ZA 0, TZA QQ, and the pressure model (21], 0, Q] We establish the 
value of scale factor at z = with a — 1. For the initial condition, we set the Gaussian 
density field with scale-free spectrum; 

V{k) oc fc n (n = -1,0,1), (11) 
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SCDM, and LCDM model. The initial condition was produced by COSMICS 

For N-body simulation, we execute P 3 M code. The parameters of simulation were given 
as follows: 

Number of particles : N = 64 3 , N = 128 3 (Figure H only) , 
Box size : L = 64/i _1 Mpc , 
Softening Length : e = 0.05/i Mpc , 
Course-Graining Length : I = 1, 2, 4/i _1 Mpc , 
Hubble patameter : h = 0.71 . 

For CDM models, we choose cosmological parameters as follows: 

SCDM : tt m = 1.0, tt A = 0.0, (7 8 = 0.84 
LCDM : tt m = 0.27, fi A = 0.73, a 8 = 0.84 

In the pressure model, we choose a polytrope exponent 7 = 4/3,5/3: In the case of 
7 = 4/3, we obtain the simplest perturbative solution given by Equation (jSJ) that is described 
by the power- law of the time variable. 7 = 5/3 is given from theoretical argumentation 
by Buchert and Dominguez Q|. They argued kinematic considerations using collisionless 
Boltzmann equation and derived 7 = 5/3. The exponent is equivalent to the adiabatic 
process of an ideal gas. Because we cannot decide on a proportional constant k in the 
equation of state from past discussion, we choose several values. In this paper, instead of k, 
we write an initial (a = 10~ 3 , i.e. z = 1000) Jeans wavenumber, given by Equation ©. 

Here we show how we set up the initial condition in the pressure model. We adjust the 
initial condition in the pressure model to be the same as that in ZA: the initial peculiar 
velocity in the pressure model is same as that in ZA. The procedure for setting up initial 



condition was shown in our previous papers 
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FIG. 1: The cross-correlation coefficient of density fields between the N-body simulation and 
Lagrangian approximations. Primordial density fluctuation is given by the scale- free spectrum 
P(k) oc k 1 (iV = 64 3 ,Z = I/i _1 Mpc). (a) When we choose 7 = 4/3, the function deviates from 
that of ZA in linear regime, (b) In the case of 7 = 5/3, we can obtain a better result than TZA. 

A. Cross-correlation coefficient 

First we calculate the cross- correlation coefficient of density fields. The Cross-correlation 



coefficient was used for the comparison of the resulting density fields 
The cross-correlation coefficient is denned by 

\ 0-10-2/ 

where <Tj means the density dispersion of model i, 



12. 
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25 



26 



27 



(12) 



(13) 



S = 1 means that the pattern of density field of two models coincide with each other. In 
linear regime, the density dispersion remains er <C 1. Although we develop the structure 
until it becomes strongly nonlinear regime (a N-body > 1); we especially analyze it in the 
quasi- nonlinear regime {a N-body — 1)- 

Figure shows a comparison of N-body density fields with those predicted by various 
Lagrangian approximations. First, we notice scale-free spectrum cases (Figure As in 
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FIG. 2: The cross-correlation coefficient of density fields between the N-body simulation and 
Lagrangian approximations (N = 64 3 ,Z = l/i _1 Mpc, scale-free spectrum model), (a) P(k) oc k°. 
The case of 7 = 4/3. (b) P(k) oc k°, the case of 7 = 5/3. (c) -P(A;) oc k . In this model, 
the difference in the coefficient between ZA and TZA becomes small. The case of 7 = 4/3. (d) 
P{k) oc fe -1 , the case of 7 = 5/3. 
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the past analyses, TZA shows better performance than ZA. Our analyses also show similar 
tendency, i.e. our analyses do not contradict past analyses. 

In the pressure model, the performance strongly depends on the polytrope exponent 7 
and Jeans wavenumber. In the case of 7 = 4/3, when we set the initial Jeans wavenumber 
to be small, even if in linear regime, the approximation deviates from an N-body simulation. 
Only for the case of Kj = 64 the approximation shows better performance than ZA in a 
quasi-nonlinear regime. We notice that the result strongly depends on the Jeans wavenumber 
in the case of 7 = 4/3: When we slightly change the value of the Jeans wavenumber, the 
cross-correlation coefficient changes dramatically. In the case of 7 = 5/3, although the result 
depends slightly on initial Jeans wavenumber, the pressure model shows a better performance 
than ZA in a quasi-nonlinear regime. Furthermore, the pressure model also shows a better 
performance than TZA. However, when we consider scale-free spectrum models, the model 
does not have typical physical scale: the model has only box size, grid size, and softening 
length. The trend of the result was unchanged when we changed the box size of the model 
and the softening parameter. When we change the number of particles, the result changes. 
From a comparison of Figure Q arid Figure El we can notice that the results depend on the 
ratio of grid size and initial Jeans wavenumber. In our calculation, we found out that it 
was good to set up the value of k so that the initial (z = 1000) Jeans wavenumber Kj was 
jV 1 / 3 /4 < Kj < iV 1//3 . For example, in the case of N = 64 3 , as we see in Figure ^ and 121 it 
is good to choose the initial Jeans wavenumber 16 < Kj < 64. 

Next we consider SCDM and LCDM models (Figure In these models, the difference 
between ZA and TZA becomes very small. Because the initial density spectrum in the CDM 
model dumps power in a small scale, the cutoff in spectrum weakly affects the formation of 
caustics, as we saw in the case of P(k) oc k~ x . From Figure EJ we can see that the effect of 
pressure improves the approximation in the quasi-nonlinear stage. Also, in both the SCDM 
and LCDM models, the case of 7 = 4/3 shows deviation from ZA in linear regime. On 
the other hand, the case of 7 = 5/3 shows that the cross- correlation coefficient becomes 
almost the same in linear regime. In the case of 7 = 4/3, when we choose a small initial 
Jeans wavenumber (for example, Kj = 16), although we can improve the approximation 
much more in quasi-nonlinear stage than in large Jeans wavenumber (Kj = 32, 64) cases, 
the approximation changes slightly for the worse in the linear stage. On the other hand, 
when we choose 7 = 5/3, although the effect seems small, we can obtain a improved solution 
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FIG. 3: The cross-correlation coefficient of density fields between N-body simulation and La- 
grangian approximations. Primordial density fluctuation is given by scale-free spectrum P(k) oc k 1 
(N = 128 s , I = l/i _1 Mpc). From comparison between Figure ^ and this graph, we can notice that 
the results depend on the ratio of grid size and initial Jeans wavenumber. 

both in linear and in quasi-nonlinear stages. 

When the model evolves to strongly nonlinear regime, the trend of the solutions change. 
In the linear stage, the case of 7 = 4/3 shows deviation from ZA. However in strongly 
nonlinear regime, though the Lagrangian approximation generally become worse, the case 
of 7 = 4/3 shows a rather good result (Figure 0] (a) and (c)). This tendency was unchanged 
even though the course-grained length was changed (Figure El and EJ) . 

In both the SCDM and LCDM models, when we choose a small initial Jeans wavenumber 
Kj, although the approximation is improved after the quasi-nonlinear stage, the reasonable 
range of Jeans wavenumbers seems wide. The strict limitation to the value of k or the initial 
Jeans wavenumber will be given by other physical considerations or by the high-resolution 
N-body simulation. 

From these results, we find that it is reasonable to choose the polytrope exponent 7 = 5/3 
until quasi-nonlinear regime is reached. These results support the suggestion by Buchert 

n 

and Dommguez |l5( and by Dommguez [19j . If we have interest in strongly nonlinear regime, 
although the Lagrangian approximation generally becomes worse, we had better analyze the 
case of 7 = 4/3. In this regime, it is necessary to consider whether that approximation will 
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FIG. 4: The cross-correlation coefficient of density fields between an N-body simulation and La- 
grangian approximations. The primordial density fluctuation is given by the CDM spectrum 
(N = 64 3 ,Z = l/i _1 Mpc). The model in which hardly a difference appears in is excluded from 
the graph, (a) The SCDM model with 7 = 4/3. (b) The SCDM model with 7 = 5/3. (c) The 
LCDM model with 7 = 4/3. (d) The LCDM model with 7 = 5/3. 
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FIG. 5: The same as Figure @J In these figures, we changed the course- grained length to I = 
2/i- 1 Mpc. (a) The SCDM model with 7 = 4/3. (b) The SCDM model with 7 = 5/3. (c) The 
LCDM model with 7 = 4/3. (d) The LCDM model with 7 = 5/3. 

still be valid. In any case, from the cross-correlation coefficient, we can give limitation in 
the polytrope exponent 7. 

Unfortunately, in these results, we cannot give a strict limit to the proportion coefficient 
k of the equation of state. When we choose 7 = 4/3, we show that the result strongly 
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FIG. 6: The same as Figure @J In these figures, we changed the course- grained length to I = 
4ft -1 Mpc. (a) The SCDM model with 7 = 4/3. (b) The SCDM model with 7 = 5/3. (c) The 
LCDM model with 7 = 4/3. (d) The LCDM model with 7 = 5/3. 

depends on k, and we can notice a strict limitation. However when we choose 7 = 5/3, we 
can hardly judge the best value for k. In our calculation, we found out that it was good to 
set up the initial (a = 10~ 3 , i.e. z = 1000) Jeans wavenumber Kj as N 1 ^ 3 /4 < Kj < N 1 ^ 3 . 
From the range of Kj, we can obtain a reasonable value for k. If we choose large value for 
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k, it becomes hard to form nonlinear structure. On the other hand, if we choose small value 
for k, the structure becomes almost the same as the structure which was obtained by ZA. 

Although the cross-correlation coefficient is one thing which is good for checking the 
accuracy of the approximation, it is not enough for checking. Now we consider two samples 
A and B. The density contrast of the samples is given by 5 a and 5b, respectively. We assume 
that the following proportion relation exists between 5a and 5b'- 

5 A (x5 B - (14) 

Even if the density contrast is greatly different, as 5a and 5b, and one shows a highly 
nonlinear structure and another remains in the linear regime, the cross-correlation coefficient 
between 5a and 5b becomes 1. 

Therefore, we must check the accuracy of the approximation with another property. In 
the next subsection, we analyze the probability distribution function of density fluctuation. 



B. Probability Distribution Function of Density Fluctuation 

Here, we compare the probability distribution function (PDF) of density fluctuation. In 
the Eulerian linear approximation, if initial data is given by a random Gaussian distribution, 
the PDF of density fluctuations will retain its Gaussianity during evolution. On the other 
hand, in the Lagrangian approximation, there appears a nonlinear effect. In fact, Kofman et 
al. |29| shows that the PDF of density fluctuation approaches a log-normal function rather 
than a Gaussian function in the cases of the Lagrangian approximation and N-body simu- 
lation. Padmanabhan and Subramanian also discussed the PDF of density fluctuation 
with the ZA and found a non-Gaussian distribution. 

How will the PDF of density fluctuation change if we take the effect of the velocity 
dispersion into consideration? Figures [71 |Hl and E] shows the PDF of density fluctuation. 
As in past works, the PDF of density fluctuation becomes log-normal in form in the N- 
body simulation. In Figures El and El the cases of 7 = 4/3 obviously show the different 
tendency: in these cases, the effect of pressure suppress the growth of positive fluctuation 
(Figure[7^b),|Hl(c),(d), andEl(c),(d)). When we also consider the PDF of density fluctuation, 
we can notice that it is not so good to choose 7 = 4/3 to examine the growth of structure, 
though the cross-correlation coefficients well show the trend. On the other hand, the case 
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FIG. 7: The PDF of density fluctuation for a scale-free spectrum (P(k) oc k, I = l/i _1 Mpc: 
^N-body - 1 at o = 1.0.). (a) The PDF of density fluctuation. In the case of 7 = 4/3, the effect 
of pressure suppresses the growth of the fluctuation, (b) The difference in the PDF of density 
fluctuation. In this figure, the difference between the case of 7 = 4/3 and other cases becomes 
clear. When we choose Kj = 32 for the case of 7 = 4/3, more greater difference appears. 

of 7 = 5/3 well show the trend in the PDF of density fluctuation. Although the difference 
of distribution between ZA and the pressure model is still small in quasi-nonlinear regime, 
the effect of the pressure can promote the evolution of nonlinear structure. Therefore the 
probability of low and high density regions increases in the case of 7 = 5/3. Furthermore, 
according to Figure Efc), the PDF of density fluctuation in the cases of 7 = 5/3 show that 
it is much better than that in the TZA case. Of course when we reach a strongly nonlinear 
regime, it is necessary to consider whether that approximation is still valid or not. 

From both the cross-correlation coefficient and PDF of density fluctuation, we can decide 
that it is reasonable to choose 7 = 5/3 as the polytrope exponent of the equation of state. 
However it is hard to decide the proportional parameter k. From the results in this paper, 
we cannot give tight limit to k. To decide the value of k, we will analyze a high-resolution 
N-body simulation or consider other physical processes. For example, we will consider the 
effect of anisotropic velocity dispersion [32] or the higher order velocity cumulant. 
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FIG. 8: The PDF of density fluctuation in the SCDM model (Z = 4/i _1 Mpc). In the case of 
7 = 4/3, the pressure effect suppresses the growth of density fluctuation. Therefore the probability 
of a small fluctuation (|<5| < 1) increases, (a) The SCDM model, at a = 0.1 (z = 9, I = 4/i _1 Mpc. 
Quasi-nonlinear regime). In the case of 7 = 4/3, the effect of pressure suppresses the growth of 
the fluctuation, (b) The SCDM model, at a = 1.0 (z = 0, / = 4/i~ 1 Mpc. Strongly nonlinear 
regime), (c) The difference in the PDF of density fluctuation between the N-body simulation and 
Lagrangian approximations. At a = 0.1. (d) Same as (c), but at a = 1.0. 
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Density distribution (LCDM, a=0.1) 
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FIG. 9: The PDF of density fluctuation in the LCDM model (I = 4/i _1 Mpc). In the case of 
7 = 4/3, the pressure effect suppresses the growth of density fluctuation. Therefore the probability 
of a small fluctuation (]6\ < 1) increases, (a) At a = 0.1 (z = 9, I = 4/i _1 Mpc. Quasi-nonlinear 
regime). The PDFs of density fluctuation seem similar to each other, (b) At a = 1.0 (z = 0, 
I = 4/i _1 Mpc. Strongly nonlinear regime), (c) The difference in the PDF of density fluctuation 
between the N-body simulation and Lagrangian approximations. At a = 0.1. (d) Same as (c), but 
at a = 1.0. 
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IV. DISCUSSION AND CONCLUDING REMARKS 



We compared two statistical quantities between an N-body simulation and Lagrangian 



approximations. In our earlier work 



22, 



231 ] . we solved the first-order perturbation equations 



in the homogeneous and isotropic background, and the second- order ones explicitly for the 
case 7 = 4/3,5/3 in Einstein-de Sitter Universe. We showed that the difference between 
the Lagrangian first-order and second-order approximations becomes small in the case of 
7 > 4/3. Therefore, in this paper we consider only the first-order perturbative solution for 
the case 7 = 4/3, 5/3. Then we carried out similar calculation with ZA and TZA to examine 
their difference from the past models. 

First, we compared these models by the cross-correlation coefficient of the density field 
between the N-body simulation and Lagrangian approximations. In scale-free spectrum 
cases, as well as in the previous analyses, TZA shows a better performance than ZA. In 
the pressure model, the performance strongly depends on polytrope exponent 7 and Jeans 
wavenumber. In the case of 7 = 4/3, when we set that initial Jeans wavenumber to be small, 
even in linear regime the approximation deviates from the N-body simulation. In the case of 
7 = 5/3, although the result slightly depends on the initial Jeans wavenumber, the pressure 
model shows a better performance than ZA in quasi-nonlinear regime. Furthermore, the 
pressure model also shows better performance than TZA. In the SCDM and LCDM models, 
the case of 7 = 4/3 shows deviation from ZA in linear regime. On the other hand, the 
case of 7 = 5/3 shows that the cross-correlation coefficient becomes almost the same in 
linear regime. When the model reaches a strongly nonlinear stage, although the Lagrangian 
approximation generally become worse, the case of 7 = 4/3 shows a rather good result. Of 
course, in this regime, it is necessary to consider whether that approximation is still valid. 

Second, we analyzed the PDF of density fluctuation. The cases of 7 = 4/3 obviously 
show a different tendency until quasi-nonlinear regime is reached: in this case, the effect 
of pressure suppresses the growth of structure. When we also consider the probability 
distribution of density, we can see that it is not so good to choose 7 = 4/3 to examine the 
growth of structure, although the cross-correlation coefficients perform well. On the other 
hand, the case 7 = 5/3 shows good tendency in the PDF of density fluctuation. Although the 
difference of the PDF of density fluctuation between ZA and the pressure model is still small 
in quasi-nonlinear regime, the effect of the pressure can promote the evolution of nonlinear 
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structure. The difference between the models of Lagrangian approximation becomes small 
when we calculate the evolution until strongly nonlinear regime is reached. From analyses 
of the cross-correlation coefficient of density field and the PDF of density fluctuation, we 
can decide that it is reasonable to choose 7 = 5/3 as the polytrope exponent of the equation 
of state. 

In this paper, we changed some values of the Jeans wavenumber Kj and undertook 
analysis. Will there be any relations between the 'nonlinear wavenumber' fc^/, in TZA and 
Kj ? The correspondence is as follows. For simplification, we consider the correspondence 
in case of a scale-free spectrum P(k) oc k n . According to the definition of a 'nonlinear 
wavenumber' in TZA, k^i is given from Equation JTUJ). In case of the scale-free spectrum 
P(k) = Ak n , the definition becomes 

1 



71+1 

From this definition, kj^L is written as 



a(t) 2 Ak%£ = 1. (15) 



k NL ~a~ 2 ^ 1+n l (16) 
For example, when we choose n = 1, k^L becomes 

k NL ~ -. (17) 

(X 

On the other hand, the Jeans wavenumber Kj in the pressure model is given from Equa- 
tion (0). When we choose 7 = 2, Kj becomes 

Kj ~ -. (18) 

a 

There are some different points to consider when we think about time evolution, although 
the relation seems to be as described above: First, in TZA, k^i affects only the initial 
spectrum. On the other hand, Kj affects the evolution of fluctuation. Second, although 
kNL obviously depends on the initial spectrum, we did not clarify the dependence on the 
initial condition of Kj. We think that a consideration of the physical process, which was 
not considered here, or the analysis of the N-body simulation, is necessary for a decision 
of Kj, i.e. k. We will have to think about the correspondence between the adhesion 
approximation and the pressure model. Buchert et al. showed how the viscosity term 
in the adhesion approximation is generated by a pressure-like force. Dommguez 17 
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discussed spatial coarse graining in a gravitating system and derived an evolution equation 
of 'adhesion approximation'. In the pressure model, we showed that the density distribution 
of the pressure model was similar to that of TZA in the previous paper 23|. The acute 
characteristic skeleton structure which appeared in the adhesion approximation could not 
be seen from the calculations in our previous paper. We will consider the relation between 
the viscosity term in coarse-grained equations and the pressure term in our model. Then we 
will analyze the correspondence between the viscosity term in the adhesion approximation 
and the proportional constant k in equation of state in pressure model. 

In this paper, we analyzed only density distribution. How will peculiar velocity distri- 
bution change with the effect of 'pressure' ? In ZA, the peculiar velocity is in proportion to 
the Lagrangian displacement. Then the growth rate of perturbation is independent of scale. 
Therefore, although the structure becomes nonlinear regime, if the initial condition is given 



29[. However, 



as Gaussian, the peculiar velocity distribution remains Gaussian all the time 
in the pressure model, the growth rate of the perturbation depends on the scale. Therefore 
the peculiar velocity distribution will deviate from Gaussian during evolution. Of course, 
the peculiar velocity distribution in an N-body simulation becomes non-Gaussian j^]]]. Does 
the effect of the pressure cause the occurrence of the non-Gaussian distribution? We think 
that the time evolution of the peculiar velocity distribution is one of the more interesting 
problems. 

In our model, we introduce the strong simplification that the velocity dispersion is approx- 
imately isotropic, i.e. the stress tensor is diagonal and a pressure-like term However, in 
general, the velocity dispersion does not remain isotropic in nonlinear regime. Until when 
is the assumption to ignore anisotropic velocity dispersion reasonable? Maartens et al. [32j 
discussed a relativistic kinetic theory generalization which also incorporates anisotropic ve- 
locity dispersion. Then they added these effects to the linear development of density inho- 
mogeneity and found exact solutions for their evolution. In a Newtonian description, though 
the equations is not generally closed, we will consider anisotropic velocity dispersion and 
the higher-order velocity cumulant and estimate their effects on the evolution of density 
inhomogeneity. 
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